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INFERRING  FINITE-TIME  PERFORMANCE 
IN  THE  M/G/l  QUEUEING  MODEL 

P.A.  Jacobs 
D.P.  Gaver 


1 .  INTRODUCTION  AND  SUMMARY 

The  M/G/l  model  for  a  single  server  approached  by  a  Poisson  (A)  stream 
of  arrivals  with  independent  but  otherwise  arbitrarily  distributed  service 
times  is  a  standard  tool  in  operations  research.  It  has  been  used  to 
model  situations  occurring  in  road  traffic,  computer  and  communication 
performance  evaluation  and  the  military,  and  has  been  modified  to 
accommodate  priority  and  server  breakdown  situations,  as  well  as  both 
heavy  and  light  traffic  situations  ( p  ~  AE  [S]  either  near,  but  below, 
unity,  and  near,  but  above,  zero). 

Although  the  M/G/l  model  is  rather  simplistic,  little  beyond  the 
moments  of  its  long-run  or  steady-state  distribution  are  easily  available 
in  anything  like  closed  (and  simple)  analytical  form  readily  evaluated 
numerically.  In  the  case  of  the  exponential  service  time  M/M/1  system, 
the  complete  analytical  solution  for  transition  probabi 1 ist ies  in  terms  of 
series  of  Bessel  functions  has  been  long  known,  and  numerical  transform 
inversion  by  Middleton  (1979)  has  provided  some  useful  tables  describing 
the  manner  in  which  steady  state  is  approached.  Recent  work  of  Abate  and 
Whitt  (1989,1988)  provides  some  handy  approximations,  particularly  in  the 
M/M/1  case.  The  work  of  Asmussen  and  Thorisson  (1988)  and  Newell  (1971, 
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1982)  furnishes  diffusion  approximations  and  heavy  traffic  results  as 
well. 

In  this  paper  we  explore  the  following  operationally  relevant 

Problem:  A  single  server  is  approached  by  a  stream  of  arrivals  to  be 
modelled  as  Poisson.  One  has  available  a  sample  of  service  times  obtained 
by  observing  the  system.  Using  the  arrival  rate  (assumed  known,  at  least 
initially)  and  the  observed  service  times,  estimate  the  expected  waiting 
time  encountered  by  a  new  arriving  customer  at  finite  time  t.  Further¬ 
more,  provide  an  assessment  of  estimation  uncertainty,  e.g.  confidence 
1 imi ts . 

The  above  is  just  one  prototypical  question  to  be  asked--it  is  the  one 
addressed  in  this  paper.  Other  questions  could  concern  the  probability 
that  an  arrival  at  t  will  wait  for  at  least  time  w,  or  that  the  customer 
waiting  time  has  never  exceeded  w  in  time  t  (starting  from  an  empty 
system) .  There  are  many  other  such  measures  of  system  performance  that 
are  relevant.  We  concentrate  here  on  inferring  the  mean  waiting  time  at 
finite  time  t,  and  do  so  non -parametrical lv ,  i.e.  without  specifying,  and 
estimating  parameters  in  a  small  conventional  parametric  family  such  as 
the  Gamma.  Our  approach  is  to  work  with  the  well-known  Laplace  trans¬ 
form  of  the  mean  (virtual)  waiting  time  for  the  M/G/l  system:  cf .  Takacs 
(1962)  and  Gaver  and  Jacobs  (1987) ,  an  empi rical  version  of  which  can  be 
obtained  by  utilizing  the  empirical  transform  of  the  busy  period,  itself 
depending  upon  the  empirical  transform  of  the  service  time  df ,  derived 
from  the  observed  data. 

We  look  upon  our  present  results  as  exploratory  and  tentative,  but  of 
an  accuracy  useful  until  refinements  are  available. 


2 


A  recent  paper  by  Tengels  (1988)  provides  mathematically  rigorous 
asymptotic  (large  sample  size)  treatment  of  certain  aspects  of  our 
nonparametr i c  infei ence  problem.  Application  of  bootstrap  methods,  cf . 
Efron  (1979)  can  yield  further  insights  into  the  inference  at  the  price  of 
computer  intensivity. 


Accession  For 


NTIS  GRA&I 
DTI C  TAB 
Unaiuicuticed 
Ju::ti  fic.itloSL 
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□  □ 


2.  MATHEMATICAL  FORMULATION 


Consider  an  M/G/l  queue  with  known  arrival  rate  A  and  independent 
identically  distributed  service  times,  S,  having  nth  moment  E  [Sn]  <  oc . 
Suppose  there  are  no  customers  either  waiting  or  being  served  at  time  0. 

Let  Vt  be  the  virtual  waiting  time  at  time  t  and  put  (conditional  on 

W0  =  0) 

f00 

t^vv(s)  =  e -st '  E  [Wt]  dt,  s  >  0  .  (1) 

J0 

We  know  that 

vw(s)  =  ~2  1  +  Poo(s)  s  ;  (2) 

where : 


Poo(s) 


_ __1_ 

S  +  A  [_1  - 


*Ks)]  ’ 


(3) 


b(s)  being  the  Laplace  transform  of  a  busy  period  duration;  b(s)  is  the 
smallest  positive  solution  of 

b(s)  =  Fs (s  +  A ( 1  -  b(s)))  :  (4) 


Fs  is  the  Laplace-St ielt jes  transform  of  the  service  time.  See  Takacs 
(1962)  and  Gaver  and  Jacobs  (1986) . 

Two  approximations  to  E  [W't  |  W0  =  0]  will  be  detailed  below;  one  for  the 
case  of  the  traffic  intensity  AE  [S]  <  1;  the  other  for  the  case  of  the 
traffic  intensity  AE  [S]  >  1.  Both  are  based  on  the  transforms  (2),  (3), 
and  (4)  . 
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3.  APPROXIMATIONS  FOR  STABLE  QUEUES :  EXPONENTIAL  APPROACH 

In  this  subsection,  it  is  assumed  that  p  =  AE  [S]  <1.  In  this  case 
there  is  a  limiting  expected  virtual  waiting  time 

^  -  frf^  •  <5) 

with 

p  =  AE[S]  ,  (6) 

the  traffic  intensity. 

Our  choice  of  the  form  of  the  approximation  is  based  on  an  observation 
of  Odoni  and  Roth  (1983),  and  on  relaxation  time  calculations  made  by 
Keilson  (1979),  Newell  (1971)  and  others.  The  explicit  approximation  of 
the  exponential: 

E;vt]  *  fx(l  -  «'*) 


or 


DC  *  E[W  *,3t 


as  e 


(7) 


It  should  be  noted  that  in  many  cases  the  exponential  approach  to  the 

limiting  value  is  not  exact,  being  only  an  approximation,  cf .  Asmussen  and 
Thorissan  (1988).  Setting 


^  _-St  E[Vt]  - 

— r: 


f~  se-st  E  dt  =  -  P  se’st  e‘5t  dt 

Jr 


,oc 

0 


(8) 


we  obtain 


st/>w(s)  -  Esq 
Coc 


— 3— 
s  +  0 


(9) 


which  is  now  solved  for  0  at  selected  s-values: 


1  +  Sii)  *  ■  s'v^TTf-  e.  ; 
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so 


a/1'1  _  ~^3C  ~  [st>w(s)  ~  ^oc  ] 
5  ~  svw(s)  -  Woe 


F i nal ly 


fi  =  -g2»w(s)  . 

svw(s)  -  Wgo 

The  approximation  to  E[Vt]  is 


E[Vt]  «  Woc(l  - 


-3t 


)  , 


(10) 


(11) 


(12) 


where  3  is  evaluated  as  above  at  s  =  ^ .  It  turns  out  that  the  ,3-values  so 
obtained  are  not  constant  but  are  (extremely)  slowly-changing  functions  of 
t,  as  should  be  the  case  from  asymptotic  analysis. 

Table  1  reports  the  values  of  the  approximating  E  [V*t]  for  an  M/M/1 
queue  with  A  =  1  and  E [S]  =  .95.  The  approximating  values  are  compared  to 
values  obtained  by  Middleton  (1979),  who  used  a  numerical  inversion  of  the 
Laplace  transform  vw  due  to  Stehfest  (1970).  The  approximation  is  very 
good  for  the  times  less  than  or  equal  to  100.  For  the  times  greater  than 
100  th**  approximation  is  within  11T(  of  the  true,  and  conservatively 
slightly  overestimate  of  the  true  mean. 
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Table  1 


Exponential  Approximation 
M/M/1  Queue:  A  =  1,  E [S]  =  .95 

T  ime 

True 

Exponent i al 

t 

E0t] 

Approx i mat  ion 

20 

3.9 

3.7 

40 

5.5 

5.4 

60 

6.6 

6.6 

80 

7.4 

7.6 

100 

8.1 

8.4 

400 

13.2 

14.6 

800 

15.5 

17.0 

1200 

16.6 

17.7 

Note  : 

=  18.05 

Table  2  reports  the  values  of  the  approximating  E[VtJ  for  an  M/G/l 
queue  with  A  =  1  and  gamma  service  times  with  shape  parameter  a  =  .2  and 
scale  parameter  3  =  4.5.  Once  again  the  approximation  is  very  good.  The 


values  of  the  true  E [Vt]  are  from  Middleton  (1979). 


Table  2 


Exponential  Approximation 
M/G/l  Queue:  A  =  1 

Service  Time  Distribution  Gamma  shape  parameter 
a  =  .2,  scale  parameter  0  =  4.5,  E [S]  =  .90 


Time 

True 

Exponent i al 

t 

E[Wt] 

Approximat ion 

4.9 

2.3 

2.2 

9. T 

3.7 

3.4 

29.2 

6.9 

6.7 

48.6 

8.8 

8.8 

97.2 

11 .9 

12.5 

194.4 

15.3 

16.9 

486 

20.1 

21  .9 

972 

22.6 

23.7 

1458 

23.5 

23.8 

1944 

23.9 

24 . 1 

Note:  =  24.3 

The  numerical  results  obtained  suggest  that  in  general  the 
approximation  is  biased  slightly  on  the  high  side  for  large  t.  It  is  thus 
conservative  in  the  sense  that  predictions  of  expected  waiting  times  tend 
to  be  slightly  overstated. 
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THE  HEAVY  TRAFFIC  EXPONENTIAL  APPROXIMATION  (HTE) 


The  most  computat ional ly- intensive  part  of  the  exponential 
approximation  is  the  numerical  evaluation  of  the  transform  of  the  busy 
period  b(s)  using  equation  (4).  This  computation  has  been  done  by  search. 
In  this  section  another  approximation  is  described  which  uses  a  heavy 
traffic  approximation  to  obtain  b(s) ,  and  is  computationally  less  arduous. 

Let  {B(t)}  be  a  Brownian  motion  with  drift  v  =  AE  [S]  -  1  and 
inf  initesir  j.1  variance  <r2  =  AE[S2]  where  S  is  a  generic  service  time.  Let 
Tx  be  the  first  passage  time  to  state  0  given  B(0)  =  x.  Put 

0(s;x)  =  E[e'sT*]  .  (13> 


where 


(15) 


Of  course  this  function  is  explicit  and  readily  evaluated,  and  no  search 
i s  needed . 

Approximating  the  virtual  waiting  time  process  by  the  Brownian  motion 
with  drift,  the  proposed  approximation  for  the  transform  of  the  busy- 
period  is 

W*)  =  E[ea(,)S]  ,  (16) 


the  Laplace-Stiel jes  transform  of  the  service  time  distribution  evaluated 
at  o(s) .  The  remainder  of  the  approximation  is  the  same  as  the 
exponential  approximation.  Table  3  shows  the  values  of  the  heavy  traffic 


approximation  for  the  M/G/l  queue  with  A  =  1  and  gamma  service  time 
distribution  with  shape  parameter  a  =  0.2  (very  long-tailed)  and  scale 
parameter  0  =  4.5.  The  heavy  traffic  approximation  is  not  as  good  as  the 
exponential  one  and  tends  to  overestimate  E  [Vt]  .  However,  the 
approximation  is  still  practically  adequate,  and  is  far  easier  to  compute 
than  is  the  “exact”  exponential  approximation  of  Section  2. 


Table  3 

Unsaturated  Heavy  Traffic  Approximation 
M/G/l  Queue:  A  =  1,  Gamma  Service  Times 
Shape  parameter  a  =  0.2,  scale  parameter  0  —  4.5,  E [S]  =  .90 


T  i  me 

True 

Heavy  Traffic 

t 

E[Vt] 

Approx i mat  ion 

4.9 

2.3 

2.5 

9.7 

3.7 

3.9 

29.2 

6.9 

7.5 

4S.6 

8.8 

9.8 

97.2 

11.9 

13.6 

194.4 

15.3 

17.9 

972 

22.6 

24.1 

1944 

23.9 

24.3 

Note:  =  24.3 
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.  THE  Bl'SYPERIQD  APPROXIMATION 


This  approach  is  based  on  approximating  the  rate  of  approach  of  the 
probability  o  an  empty  system,  PooC^)*  ‘to  its  eventual  value.  Since  the 
process  of  alternating  busy  and  idel  periods  can  be  expressed  in  terms  of 
an  alternating  renewal  process,  methods  of  Gaver  and  Jacobs  (198S)  can  be 
applied. 

Note  : 


s  •  ^w(s)  = 


+  PooO)  » 


where 


Poo(s)  =  I 


1  ,  , 


Poo(»)  =  f  «•**  P{x(t)  =  0|X(0)  =  0}d1; 

J  n 


where  X(t)  is  the  number  of  customers  waiting  or  being  served  at  time  t. 
Assume  the  queue  is  stable;  then 

p<x<t>  -  0|X(0>  -  °>  ■  E[IduPrlWl  =  1  -  -  •  <2°! 

where  Idle  is  the  length  of  an  idle  period  and  Busy  is  the  length  of  a 

busy  period.  Approximate  as  follows  focussing  on  the  probability  that  the 
system  is  empty,  approximates  as  follows  in  terms  of  the  parameter  0b: 

P{X(t)  =  0|X(0)  =  0}  *  (1  -  p)  +  pe'0bt  ;  (21) 


from  this  comes 


E[Vt]  *  (p  -  l)t  +  (1  - 


P) t  +  f  (1  -  e'  b  ) 


=  f  U  - 
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To  find  /?b : 


_  P 

PooCs)  =  - s -  +  j — 


T  S 


pm(s>  =  1  i ;  =  J  i  +  F±^> 


where 

c<»>  -  [1  .'em*'3  *"d  >£W  =  r~h 

Solving  for  0b  results  in 


(23) 

(24) 


(25) 


sC(s)  [j — ^ — p  +  1]  sC(s)[r4-^] 

1  -  C(s)  “  1  -  C (s ) 


(26) 


and  evaluate  at  s  =  A. 

Table  4  shows  values  for  the  busy  period  approximation  for  an  M/M/1 
queue  with  A  =  1  and  p  =  .95.  The  approximate  values  are  below  the  true 
values  for  most  times  t.  This  approximation  is  not  as  accurate  as  is  the 
exponential  approximation  of  Section  2,  no  longer  being  conservative. 

Table  4 

Busy  Period  Approximation 
M/M/1  Queue 


A  -  1,  E[S]  =  .95 

T  ime 

T  rue 

Approximate 

t 

E  [Vt] 

E[Vt] 

20 

3.9 

4.1 

40 

5.5 

5.4 

60 

6.6 

6.3 

80 

7.4 

7.0 

100 

8.1 

7.5 

400 

13.2 

11  .6 

800 

15.5 

13.5 

1200 

16.6 

14.6 
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6 .  THE  UNSTABLE  OR  SATURATED  QUEUE  EXPONENTIAL  APPROX I MAT I ON 


Rewriting,  we  have 


V'w  (  s ) 


“5 


1  + 


A[l_^b(s)] 


(27) 


=  P~7rl  +  *  s  ~  ACl1-  b(i)]  *  (2*) 

In  this  subsection  we  assume  that  p  >  1.  When  this  is  true  the 
distribution  of  the  length  of  the  busy  period  is  not  honest,  and 
b (0)  <  1.  We  will  write  b(0)  b#(s)  =  b(0)  [g^]  .  where  b#(s)  is  now  the 

transform  of  an  honest  random  variable.  Now 


w  (  s  ) 


P  ~  1  +  1  _ 1 _ 

*S  g  s  +  A  [1  -  b(0)  b*(s)] 


(29) 


For  smal  1  s 


*'w(s)  *  e-jJ  *  i  A  [l  -’btoj]  •  (3°) 

Thus  for  large  t 

E[Wt]  *  (p  -  l)t  +  A-r!  -^(0)]  •  (31) 


Further 


vw  (  s  )  ~  s2 —  +  s  ^  (  s  ) 

where 


(32) 


^(S)  =  - i -  ~  _ 1 _  , 

s  +  A[1  -  b(0)b#(s)]  s  +  A  -  Ab(0)b#(s) 


Multiplying  through  by  s  +  A  yields 


(•  +  A)*(s) 
or 


_ 1 _  . 

1  -  j-i-jb(0)b*(.)  ’ 
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(s  +  A)d(s)  =  1  +  (S  +  A)o(s)[r^-1;b(0)b#(s)]  ;  (33) 

and 

«(s)  =  "I  A  +  s~~+ — A  b(°)b#(s)^(®)  •  (34) 

r°° 

Letting  ^(s)  =  e'$t  f(t)dt 
J0 

X  f°° 

f(t)  =  e"At  +  f  (t-u) h (u)du 
J0 

where 

7(s)  =  J°°  e*sth  (t )  dt  =  b  (0)  j-A-j.  b*(s)  (35) 

and 

h (0)  =  b(0)  <  1  . 

Thus  f  satisfies  a  renewal-type  equation  with  a  degenerate  distribution 
and  it  can  be  shown  (cf.  Feller  (1966)) 

-3\ 

f  (t)  s:  ce  as  t  —  oc  .  (36) 

Thus .  since 

VW  ( s )  =  —  +  5  ^  ( s ) 

E[wtl  **  (P  ~  l)t  +  |(1  -  e-^1)  (37) 

as  t  —  oc .  Comparing  (36)  with  (31)  indicates  that 

f  =  A^l  -1 b(0)]  -  (38) 
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In  the  saturated  case  the  exponential  approximation  is  computed  as 
fol lows : 


E[Vt]  -  (p  -  i)t  »  a  (1  -  e_/Jt) 


where 


a 


*  ATT 


-  b(0)j 


Thus 


E [Wt]  -  (p  -  l)t  -  a  -0t 

g-  —  _e 


(39) 


Taking  Laplace  transforms 

s%(s)  -  (-  i--)  -  a 
- s - 


1 

1  +  /?(i) 


(40) 


which  leads  to 


-s2ww(s)  +  (p  -  1) 
svw(s)  -  (?■  g . 1)  -  1 


(41) 


s  i  nee 


a 


A [1  -  b(0)]  • 


(42) 
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THE  SATURATED  QUEUE  HEAVY  TRAFFIC  EXPONENTIAL  APPROX I MAT I ON 


Once  again  the  most  arduous  computation  in  the  saturated  queue 
exponential  approximation  is  the  computation  of  b(s),  the  transform  of  the 
length  of  the  busy  period.  As  before,  an  approximation  is  proposed  which 
uses  a  heavy  traffic  approximation  to  evaluate  b(s)  and  b(0) . 

Specifically,  since  u  =  AE[S]  -  1  >  0  let 

®(0)  =  =3?  ‘  (43) 

a 


and 


a(s) 


- V 


2<t2s 


(44) 


and  put 


b  (0)  =  E[ea(o)s] 


(45) 


and 

b(s)  =  E[eQ($)s]  . 


(46) 


The  remainder  of  the  approximation  is  as  in  the  exponential  approximation 
for  the  saturated  queue. 

Table  5  shows  results  for  both  the  “exact”  exponential  approximation 
and  heavy  traffic  approximation  for  a  saturated  M/G/l  queue  with  A  =  1  and 
gamma  service  times  having  shape  parameter  0.2  and  scale  parameter  7.5  and 
E[S]  =  1.5.  The  true  values  are  from  Middleton  (1979).  As  is  expected 
the  exponential  approximation  gives  values  which  are  closer  to  the  true. 
Both  approximations  improve  as  t  increases.  The  exponential  approximation 
is  very  good  for  all  times  except  the  very  smallest. 
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Table  5 


Saturated  M/G/l  Queue 
A  =  1,  Gamma  Service  Times,  E [S]  =  1.5 
Shape  parameter  0.2,  Scale  parameter  7.5 


T  ime 

True 

Exponential 

Heavy  Traffic 

t 

E[Vt] 

Approximation 

Approximation 

1.6 

3.9 

1.9 

2.1 

8.1 

7.5 

7.4 

8.2 

16.2 

12.9 

13.1 

14.4 

27.0 

19.3 

19.7 

21 .6 

43.2 

28.3 

28.8 

31  .3 

54.0 

34.0 

34.6 

37.3 

108.0 

61 .9 

62.2 

65.6 

162.0 

89.1 

89.2 

92.8 

216.0 

116.2 

116.2 

119.8 

324.0 

170.2 

170.2 

173.8 

432 . 0 

224.2 

224.2 

227.8 
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8.  INFERENCE 


The  approximations  to  E[Vt]  can  be  used  for  inferential  purposes  by 
replacing  moment  and  transforms  by  sample  moments  and  empirical 
transforms.  More  specifically,  suppose  X  =  1  is  known  and  service  time 
data  dlt...,dn  are  collected.  The  empirical  Laplace  transform  of  the 
service  time  distribution  is 

*s(s)  =  £  e  ld|  ;  •  (47) 

i=l 

Woe  is  estimated  by 


^  oc  — 


A? 


2(1  -  p) 


(48) 


where 


d  - 


(49) 


d3  = 


(50) 


and 


P  =  A3  .  (51) 

I 

The  empirical  Laplace  transform  of  the  service  times  and  the  moment 
estimator  of  p,  p,  can  be  used  in  formulas  (2)-(4)  to  obtain  an  estimate 
of  iw(s)'  This  estimate,  V>w(s),  together  with  W  can  be  used  in  the 
approximations  to  obtain  estimates  of  E[W(t)] . 

Tables  6-8  report  results  of  simulation  experiments  to  study  the 
behavior  of  the  estimates.  Each  simulation  has  1000  replications.  Each 
replication  supposes  a  sample  of  100  service  times.  The  random  numbers 
were  generated  using  LLRAND0MII  random  number  package,  see  Lewis  and  Uribe 
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Table  6 


Moments  of  Estimates  of  E[Vt] 
for 

M/M/1  Queue  with 
A  =  1  E[S]  =  0.9 


T  ime 

True 

Exponential 

Heavy 

Traf f i c 

t 

E[Vt] 

Approximation 

Approximation 

Mean 

Variance 

Mean 

Variance 

6.5 

1.9 

1.8 

0.1 

2.0 

0.1 

32.4 

4.2 

4.3 

2.0 

4.5 

2.4 

97.2 

5.9 

6.8 

14.0 

7.1 

15.6 

162.0 

6 . 7 

8.3 

33.1 

8.6 

35.9 

486.0 

7.9 

13.0 

210.6 

13.2 

218.5 

64S  .0 

S.O 

14.7 

347.6 

14.9 

357 . 5 

Table 

T 

Moments  of  Estimates  of  E[Vt] 
M/G/l  Queue 

A  =  1  ECS]  =  0.5 

Gamma  Service  Times 

Shape  =  0.2  Scale  =  2.5 

T  ime 

True  Exponential 

Heavy 

Traf f i c 

t 

E[Wt]  Approximation 

Apprcximat ion 

Mean 

Variance 

Mean 

Variance 

0.3 

0.13  0.12 

0.001 

0.13 

.001 

1.8 

0.53  0.51 

0.024 

0.56 

.03 

3.0 

0.71  0.70 

0.058 

0.78 

.08 

4.8 

0.90  0.91 

0.13 

1.0 

.  17 

9.0 

1.2  1.2 

0.33 

1  .3 

.44 

24.0 

1.4  1.6 

1.0 

1.6 

1.2 
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Table  8 


Moments  of  Estimates  of  E[Wt] 
M/G/l  Queue 
A  =  1,  E[S]  =  1.5 
Gamma  Service  Times 
Shape  =  0.2  Scale  =  7.5 


Time 

True 

Exponential 

Heavy 

T  raff i c 

t 

E[Wt] 

Approximation 

Approximation 

Mean 

Var i ance 

Mean 

Variance 

1.6 

3.9' 

1.9 

0.30 

2.0 

.31 

8.1 

7.5 

7.3 

6.7 

8.1 

7.6 

54.0 

34.0 

34.5 

256.8 

36.8 

277.7 

162.0 

89.1 

90.1 

2436.0 

92.9 

2488.0 

324.0 

170.2 

171.3 

10,210 

174.3 

10,296 

(1981).  If  the  true  average  service  time  is  close  to  1,  then  it  is 
clearly  possible  for  the  sample  traffic  intensity  p  to  be  less  than  or 
greater  than  1.  If  p  is  less  than  1,  then  the  estimate  using  the 
approximation  for  the  stable  queue  is  computed.  If  p  >  1,  then  the 
estimate  using  the  approximation  for  the  unstable  queue  is  computed.  This 
choice  appears  to  be  natural  unless  other  information  is  available,  or 
more  assumptions  are  made. 

In  Tables  6-8  are  reported  means  and  variances  of  the  estimates  of 
E  [Vt  ]  for  the  estimates  based  on  the  exponential  approximation  and  the 
heavy  traffic  approximation. 

In  Table  6  results  for  a  M/M/1  queue  with  E [S]  =  .90  are  given.  The 
same  random  numbers  were  used  to  compute  the  estimates  for  each  of  the 
times  t.  Of  the  1000  replications,  137  had  p  >  1  and  so  the  unstable 
queueing  approximations  were  used  in  these  cases.  Doubtless  it  is  the 
contribution  of  these  cases  that  lead  to  the  pronounced  over  estimate  of 
the  mean . 
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In  Table  7  are  reported  results  for  an  M/G/l  queue  with  A  =  1  and 
gamma  service  times  with  E  [S]  =  .5  having  shape  parameter  .2  and  scale 
parameter  2.5.  None  of  the  1000  replications  has  p  >  1. 

In  Table  8  are  reported  moments  for  the  estimates  of  E[W{]  for  an 
unstable  M/G/l  queue  with  A  =  1  and  gamma  service  times  with  shape 
parameter  0.2  and  scale  parameter  7.5  and  E [S]  =  1.5.  Of  the  1000 
replications  p  >  1  for  938.  of  them.  The  true  values  of  E[Vt]  in  all  the 
tables  are  from  Middleton  (1979) . 

The  means  of  the  heavy  traffic  approximation  are  larger  than  those  for 
the  exponential  approximation.  The  variances  increase  as  t  increases. 

The  very  large  variances  in  Tables  6  and  8  are  attributable  to  those 
estimates  of  E  [Vt]  for  which  p  >  1. 

In  order  to  put  the  large  variances  of  Table  8  into  perspective . 
consider  the  following  approximation. 

As  a  first  approximation,  if  AE  [S]  >  1,  then  E[Vt]  ss  (AE[S]  -  l)t 
which  caa  be  estimated  using  the  mean  of  the  service  times 

m(t)  =  (A3  -  l)t  . 

The  variance  of  m(t)  for  a  sample  size  of  100  is 
Var  [m  (t )  ]  =  A2  V  *  qqS  ^  t2  . 

In  the  case  of  a  gamma  service  time  distribution  with  shape  parameter  0.2 
and  scale  parameter  7.5,  A  =  1,  t  =  162 

Var [m (t) ]  =  2952.5  . 

Comparing  this  number  to  the  corresponding  variance  of  the  estimate  in 
Table  8  of  2436  indicates  that  the  latter  variance  is  not  unreasonable. 
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The  means  of  the  estimates  in  Table  7  are  close  to  the  theoretical 
values  in  this  stable  queue.  The  means  of  the  estimates  in  Table  6  are 
close  to  the  theoretical  values  for  the  smaller  t’s.  For  larger  t  the 
means  are  greater  than  the  theoretical  values  owing  to  those  replications 
for  which  p  >  1 . 

The  means  of  the  estimates  in  the  unstable  queue  case  of  Table  8  are 
close  to  the  true  values  of  E[Wt]  except  for  the  smallest  time,  t  =  1.6. 
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9.  CONCLUSIONS 


This  paper  proposes  an  easily-computed  approximation  to  the  finite¬ 
time  expected  waiting  time  for  an  M/G/l  system  starting  from  an  empty 
condition.  Both  unsaturated  (p  <  1)  and  saturated  (p  >  1)  conditions  are 
considered.  Numerical  evidence  is  presented  to  indicate  that  the  quality 
of  the  approximation  is  usefully  good,  especially  when  ease  of  computation 
is  an  issue.  Further,  the.  methodology  is  adapted  to  assess  expected 
waiting  time  when  inferences  must  be  made  from  a  random  sample  of  service 
times,  and  the  decision  is  made  to  do  so  nonparametr ical ly ,  i.e.  without 
fitting  a  specific  function.  The  results  appear  reasonable  and 
potentially  useful,  and  are  not  burdensome  to  obtain.  The  methodology 
investigated  can  also  be  applied  to  the  variety  of  queueing  models  that 
are  close  siblings  of  M/G/l:  priority  and  breakdowns  and  ’’vacations’’ 
being  examples.  Of  course  other  approximating  and  inferential  options 
remain  to  be  investigated. 
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